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ABSTRACT 

We show that outside-in disc instabilities, which can produce rapid-rise FU Orionis 
outbursts in T Tauri systems, are a natural consequence of the existence of protoplan- 
etary/protostellar companions. The protoplanetary/protostellar companion would be 
formed through gravitational instability on a dynamical timescale (~ 10 5 y), and pro- 
vided it was formed at less than a few AU, it would be swept in to ~ 15i?Q in less 
than 10 5 y. A companion of mass 10 _2 Mq at this radius would act as a flood gate 
which stores material upstream and periodically releases it following the ignition of 
the thermal ionisation instability. 
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£2 1 INTRODUCTION 

o 

I I Models of thermally unstable accretion discs have been suc- 
cessful in explaining the chief spectrophotometric properties 
Q"^°f FU Orionis outbursts in young stars ( Bell etal. 1995, Bell 
Lin 1994, Clarke, Lin & Pringle 1990). In such models, 
i~H the inner disc undergoes a thermal runaway, associated with 
*jp"the partial ionisation of hydrogen, once its surface density 
^achieves a critical value, Ea- This thermal runaway heats 
V-j the inner disc (out to a few tenths of an AU) and boosts 
X/l the local accretion rate by about an order of magnitude rel- 
ative to its value in the outer disc, from ~ 10 5 Mqj J to 
~ 10 -4 MQy -1 . In this state, the modelled disc spectral 
energy distribution is a good fit to that observed in FU Ori- 
onis systems. The models also predict the observed line 
width-excitation relation (Hartmann & Kenyon 1985). De- 
cline from outburst results from a downward thermal run- 
away that ensues once the surface density in the inner disc 
falls below another critical value. 

Reproducing the observed light curves in FU Ori ob- 
jects with such models has been more problematical how- 
ever. In a disc that is isolated from its surroundings, apart 
from mass input at its outer edge, thermal instability first 
occurs at the inner edge and then propagates outwards. The 
resulting light curve rises slowly (on a timescale of a number 
of years) and in this respect resembles the slow rise system 
V1515 Cyg. In order to achieve a rapid rise (on a timescale 
of less than a year), as in V1057 Cyg and FU Orionis, it is 
necessary that the outburst is initiated at some larger radius 
(Bell etal. 1995, Clarke, Lin & Pringle 1990), and that the 



wave of runaway heating propagates inwards as well as out- 
wards. The different temporal behaviour in this case results 
from the different ease with which such heating waves prop- 
agate inwards and outwards (see, for example, discussion in 
Lin, Papaloizou & Faulkner 1985). The triggering of the in- 
stability at larger radius can be achieved by the application 
of ad hoc density perturbations to the disc (Bell etal. 1995, 
Clarke, Lin & Pringle 1990). Whereas the resulting light 
curve then provides a very satisfactory fit to rapidly rising 
systems, the ad hoc nature of the perturbations is clearly 
undesirable, especially as no obvious mechanism has so far 
presented itself (see discussion in Section 4). 

In this paper, we point out that the triggering of out- 
side in outbursts (as required in rapid-rise systems) is a nat- 
ural consequence of the existence of a close protoplanetary/ 
protostellar companion. Such a satellite can carve out an 
annular gap around itself, which separates material in the 
inner disc from that upstream of the satellite. The satellite, 
if sufficiently massive, causes material to accumulate up- 
stream, whilst allowing the inner disc to drain away onto the 
central star. As a result, thermal runaway is first triggered 
just upstream of the satellite. As the disc heats up during 
the thermal runaway, it may be able to overwhelm the gap 
around the satellite and to flow into the inner disc. From 
this point onwards, the outburst would evolve as though 
the satellite were not there. After the outburst, the satellite 
would be able to re-create a gap in the cool disc material and 
the cycle would repeat itself, the satellite again controlling 
the location at which outburst is initiated. 

In the scenario outlined above, the satellite behaves like 



2 C.J. Clarke & D. Syer 



a flood gate, which is closed during the inter-outburst pe- 
riod and which is suddenly opened, releasing the stored up 
material, when thermal instability occurs in the material 
upstream. In what follows, we quantify this idea, making 
use of the calculations of coupled disc-satellite interactions 
of Syer & Clarke (1995). These calculations are applicable 
only to rather low mass satellites (M s < 1O" 2 M ). Calcu- 
lations in the case of more massive satellites demonstrate a 
qualitatively similar picture, in the sense that the satellite 
again generates a gap around itself that separates the inner 
from the outer disc (Artymowicz et al 1991, Artymowicz & 
Lubow 1994). The secular disc response in this case has 
been calculated only for the case that the disc is not subject 
to an external mass input, and also that the star-satellite 
system is not allowed to evolve in response to interaction 
with the disc (Pringle 1992). We stress that we concentrate 
on the case of lower mass satellites purely because the neces- 
sary calculations are available, but that we would expect the 
idea to work in a qualitatively similar manner for a binary 
companion more nearly equal to the primary mass. 

The structure of the paper is as follows. In Section 2 
we describe the coupled evolution of a satellite-disc system 
and calculate the radius at which thermal runaway is first 
triggered, comparing this radius with the requirements of 
time-dependent disc models for FU Orionis outbursts. In 
Section 3 we discuss the nature of the companion, and the 
constraints that are placed on its formation and evolution. 
Section 4 is a discussion. 



2 DISC-SATELLITE EVOLUTION. 

A satellite, mass M s , orbiting a star of mass M„, clears out 
material in an annular gap whose extent is set by a balance 
between tidal and viscous effects. The gap remains open 
only if it is larger than both the satellite's Roche radius 
and the vertical pressure scale height, H, of the disc (Lin & 
Papaloizou 1985). In the case of viscous a discs (Shakura 
& Sunyaev 1973) these criteria demand that the disc aspect 
ratio (H/r, where r is the distance from the central star) at 
the planet should satisfy: 



H/r < min 



1/2 



1/5' 



(1) 



Syer & Clarke (1995) computed this coupled satellite- 
disc evolution in the limit of a planet relatively massive com- 
pared with the local disc, and showed that the results could 
be understood in terms of a simple analytical model. In this 
model, the ratio of the surface density just upstream of the 
satellite, E s , relative to its unperturbed value at that point, 
Eo, depends on a single parameter B: 



B 



47rS r s 2 



(2) 



\Beq 



where r s is the instantaneous radius of the satellite. B is 
thus a measure of the satellite mass relative to the local 
unperturbed disc mass: the analysis presented here is appli- 
cable in the limit B < 1. If the relationship between Eo and 
the accretion rate, M, in a steady state disc is given by the 
power law: 



Eo oc M a r b 



(3) \Soprop 



then the surface density enhancement at the satellite is given 
by: 



E s = E 5 



-a/(l + a) 



( 4 ) \Sseq 

and the migration timescale of the satellite is related to the 
unperturbed local viscous timescale by: 



is = toB 



-a/(l + a) 



(5) 



\tmeq 



(Syer & Clarke 1995) where q = M s /M t and a is the 
usual viscosity parameter. Once a satellite has opened up a 
gap, its evolution is coupled to the disc by differential tidal 
torques across the gap. Depletion of the inner disc by ac- 
cretion onto the star results in an unbalanced tidal torque 
on the satellite, which therefore spirals inwards. In the limit 
that the satellite is light compared with the local disc (in the 
sense defined below) the satellite behaves like a representa- 
tive fluid element of the disc, and therefore migrates inwards 
on the disc's viscous timescale. A heavier satellite is unable 
to migrate on this timescale because of its higher inertia. 
As a result, material accumulates behind the satellite until 
it can extract enough angular momentum from the satellite 
to drive it inwards. The disc surface density profile adjusts 
itself so that its edge moves inwards with the satellite and 
so that angular momentum is conserved globally. 



(Syer and Clarke 1995). 

The nature of the tidal coupling between disc and satel- 
lite falls into various regimes, separated by three critical 
radii: r is the radius at which the gap closure criterion 
(equation 1) is marginally satisfied; tb, is the radius at 
which 5 = 1 (equation 2); and ta is the radius at which 
E s = Ea, and thermal runaway ensues. Since H/r in a 
steady state disc is generally an increasing function of r, it 
follows that the satellite can open up a gap when r s < r . 
Inward of r , the satellite is tidally coupled to the disc and 
spirals inwards on the local viscous timescale. For r s < tb, 
the disc behind the satellite (r > r s ) is dammed up according 
to equation (4). The downstream disc (r < r s ) can accrete 
onto the central star, so that a disparity develops between 
the surface densities inward and outward of the satellite. If 
XHrmi^g va j ue Q f ff j r behind the satellite at ta is boosted so as to 
exceed the limit given by equation (1), the gap around the 
satellite is overwhelmed and heating waves propagate both 
in and out. Provided that ta is significantly larger than 
the disc inner edge, the character of the ensuing outburst is 
'outside-in' (and therefore rapidly rising, Bell et al 1995). 

We now evaluate tb and ta for plausible models for 
the disc structure around young stellar objects. For the low 
values of a invoked in such discs (10 -4 — 10 -3 , Bell & Lin 
1994) and for M s of a giant planet mass or above, the critical 
value of H/r for gap closure is large, in the range several 
tenths to 1. Since the disc is thin between outbursts (H/r < 
0.1, Bell et al 1995), it follows that such a satellite can open 
a gap at all radii. In order to calculate tb, it is necessary 
to prescribe the radial dependence of the surface density 
in the unperturbed disc, which we parameterise (from the 
equilibrium solutions of Bell & Lin 1994) as the power law 
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Figure 1. The surface density of the disc £ as a function of 
r (in solar radii). The dot-dashed line is the unperturbed disc, 
So . The long-dashed line shows the critical surface density Y^a 
for triggering the ionisation instability. The solid line shows the 
effect of a satellite in banking up the disc S s (r s ) . For parameters 
see text. 



E =6 x 10 4 gcm" 2 

/ r x -Q 8 / a x -Q 8 / M \° 65 ( 6 ) 

vau; uo- 3 ; ^lo-^Qy- 1 J 

where we have used a = 0.65 and b = —0.8 (equation 
and thus 



10 -6 Mqv -1 and a = 10 -3 . The solid line depicts the bank- 
ing up of the surface density inward of tb according to equa- 
tion (4), whilst the dashed line represents the radial depen- 
dence of the critical surface density Ea (equation (8)). We 
find that the location of ta depends rather weakly on M , a 
or M s and lies in the range .02-.1AU (5-25_Rq) for a satellite 
in the mass range 10 - -10 - Mq. A satellite of Jupiter mass 
or above would bank up the disc behind it so that thermal 
instability would be triggered at a radius of around 10i?Q. 
The value of H/r once the instability is triggered is large 
(several tenths, Bell & Lin 1994), so that one would expect 
the gap to be overwhelmed once the thermal runaway is un- 
der way. The instability would then sweep into the inner 
disc and produce a rapid rise FU Orionis outburst. 

The outburst is terminated when the disc surface den- 
sity falls below a critical value, triggering downward thermal 
runaway to a cool state. At this point, the disc is thin again, 
so that the satellite can re-open a gap. The surface density 
in the post-outburst disc is low, due to draining of mate- 
rial onto the star during outburst, so that the effect on the 
satellite orbit is small at first. In fact, the surface density 
\gap required to resume pushing the satellite downstream is given 
by equation (4) and is therefore, by construction, equal to 
Ea at ta- Thus the satellite, migrating inwards through 
tidal coupling to a disc, is 'stalled' at the radius at which 
thermal instability is triggered. The reason for the stalling 
is that at this radius material periodically overflows past 
the satellite and therefore relieves the torque that otherwise 
would drive it inwards. We therefore have the possibility 

. „ that a satellite can remain at ta over a substantial period, 
\Sono . . . . 

acting as a gate for generating repeated FU Orionis out- 
bursts. The constraints that this requirement places on the 
satellite are discussed in Section 3 below. 



r B =0.4AU 



io- 2 M e J Vio 



M 



(7). 



Adopting the following radial dependence for the critical 
surface density for the triggering of the thermal ionisation 
instability (again from Bell & Lin 1994): 



Ea = 10 g cm 



(JL-Y 



(8) 



vau; vio- 3 

we obtain the result that the surface density at the satellite 



E s =3 x 10 gem 



4 -2 



/ M \° V M s (9) ' 



.10- 3 ; ^lO^Mgy" 1 J \io- 2 M e 
(equations (4) and (6)) equals Ea at a radius 



r A =0.06AU 



10- 2 M, 



© 



10- 



M 



lO^Mgy- 1 



(10) 



Figure 1 illustrates the surface density profile of a 
steady state disc (dash-dotted line, equation 6) with M = 



3 FORMATION AND EVOLUTION OF THE 
SATELLITE 

\rBnoFU Orionis outbursts are believed, on statistical grounds, to 
occur ~ 100 times over the lifetime of a T Tauri star (if all 
T Tauri stars undergo FU Orionis outbursts). Furthermore, 
the phenomenon of multiple bow shocks in Herbig- Haro 
objects has been tentativley interpreted as resulting from 
FU Orionis outbursts at intervals of ~ 1000 years. This 
leads to a period over which FU Orionis outbursts occur of 
\SAnoa few times 10 5 years. Since the lifetime of the T Tauri 
phase is much longer than this, it is natural to associate 
FU Orionis outbursts with the early epoch of the T Tauri 
phase. This association is consistent with the high mass- 
input rates (~ 10 -5 MQy -1 ) required in disc models of the 
outbursts (Kawazoe & Mineshige 1993, Bell & Lin 1994) 
and also with the observed nebulosity, indicative of youth, 
around bursting objects (Kenyon et al 1987). If FU Orionis 
outbursts are indeed confined to young systems, then the 
satellites that we are postulating must have formed quickly. 

This requirement presents no problems in the case of 
satellites formed through gravitational instability, either on 
observational or theoretical grounds. Binary systems are 
observed amongst the youngest T Tauri stars (age ~ 10 5 y, 
Mathieu 1994), indicating that companions can form on a 
timescale comparable with the free-fall time of the parent 
molecular cloud core. Extensive numerical simulations of 
collapsing rotating clouds also demonstrate the possibility 



\ssno 
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of binary fragmentation on a dynamical timescale (e.g. Boss 
& Bodenheimer 1978, Boss 1987, Bonnell et al 1991, Bate 
etal. 1995). The low mass companions that we have consid- 
ered in our calculations in Section 2 are intermediate in mass 
between the giant planets and stellar objects, and are too 
small for their formation to be resolved in the fragmentation 
calculations mentioned above. Such rapid formation is more 
problematical if the satellite is assumed to form by prior ag- 
glomeration of a rocky core, where timescales of ~ 10 8 years 
are usually considered necessary (Bodenheimer & Pollack 
1986, Nakagawa et al. 1983). This latter estimate is how- 
ever inconsistent with the low gas masses observed around 
pre-main sequence stars (Zuckerman et al 1995) which im- 
plies that rapid formation (as in the models of Lissauer 1987 
and Barge & Sommeria 1995) is required. We stress that al- 
though the satellites we postulate overlap the giant planets 
in mass, they do not necessarily resemble them in structure 
and formation history: the satellites considered here are ul- 
timately accreted by the central star and so do not survive 
the pre-main sequence stage. 

The second requirement of our model is that the com- 
panions, once formed, can find their way to a radius ta 
(~ 0.1AU) within ~ 10 5 years. This requirement is sat- 
isfied for objects that form at sufficiently small radii that 
disc-satellite coupling can bring them in on the required 
timescale. For r s > tb (equation 7) the satellite spirals 
in on the disc's viscous timescale 

t =4 x 10 4 y 

/ r \ 12 / a \ 08 / M V 35 ( n ) 
VAU7 V10- 3 J ^lO^Mgy- 1 J 

as derived from the equilibrium solution (6). For r s < 
tb the migration timescale is given by equation (5): al- 
though this implies that t s increases relative to to at small 
radii, the steep inward decrease of to means that the largest 
timescales are still encountered at large radii. Thus, from 
equation (11), the requirement that the satellite can be 
swept into a radius ta within ~ 10 5 years translates roughly 
into a requirement that it be formed at radii of a few AU or 
less. 

The final requirement of our model is that the satellite 
is stable against tidal disruption at r s = ta, which effec- 
tively means that it must already have undergone the hydro- 
dynamic collapse associated with dissociation of molecular 
hydrogen at its core. This point is reached on a Kelvin- 
Helmholtz timescale when the central temperature reaches 
2000 K. (Note that a satellite that is disrupted can only 
give rise to a single FU Orionis outburst, caused by its mass 
accreting onto the central object. Whereas one might pos- 
tulate a wealth of protoplanetary/protostellar companions 
so as to explain FU Orionis outbursts as multiple satellite 
disruptions, we concentrate here on the more economical 
solution of using a single satellite). 

From models of the formation of brown dwarfs (Burrows 
etal. 1993), we find that an object of mass ~ O.OIMq at an 
age of around 10 5 years is more dense, on average, than a 
T Tauri star. Such an object would fill its Roche lobe only 
if brought close to the surface of the T Tauri star, and thus 
is not susceptible to tidal disruption at r s = ta- A simi- 
lar conclusion can be reached from models of the formation 



of Jupiter through gravitational collapse (Bodenheimer et 
al 1980). Note that before hydrodynamic collapse the pro- 
toplanetary/protostellar companion must be stable against 
tidal disruption at its site of formation. This condition is sat- 
isfied in models of formation of giant planets through grav- 
itational instability at about 5AU(e.g. Bodenheimer etal. 
1980). 

Following hydrodynamic collapse, which occurs on a 
timescale of around a year, the resulting structure is com- 
pact, with a high central temperature (~ 30, 000 K). It is 
notable that the high central temperature of protoplane- 
tary/protostellar companions would protect them from com- 
plete evaporation during the FU Orionis outburst, when im- 
mersed for periods of ~ 100 years in gas at several times 
10 4 K. Partial evaporation of the outer layers would be ex- 
pected however: the temperature inversion induced by heat- 
ing the satellite from the outside would tend to supress con- 
vective cooling and thus puff up the outer layers, as energy 
released from the satellite's gravitational contraction would 
be unable to escape (Cameron etal. 1982). Substantial (i.e. 
order unity) changes in the satellite's radius would how- 
ever requires that such immersion lasted a Kelvin-Helmholtz 
timescale. Thus the brief interludes of immersion envisaged 
here would be expected to have only a modest influence on 
the tidal stability criterion of the satellite. 



4 CONCLUSIONS 

We argue that, in the scenario outlined above, the presence 
of the satellite sets up the disc in a state suitable for the 
generation of outside-in outbursts, as required in rapid rise 
FU Orionis systems. This argument is based on the detailed 
investigation of Bell etal. (1995) who found that essentially 
any perturbation that destabilised the disc at radii signif- 
icantly greater than its inner edge gave rise to outside-in 
outbursts of the required form. Past conjectures about the 
origin of such perturbations have involved the effects of a 
massive binary companion on an eccentric orbit that peri- 
odically disturbs the disc (e.g. Kenyon et al 1988, Clarke et 
al 1990, Bonnell & Bastien 1992, Bell etal. 1995), but the 
absence of observed companions at suitable radii is a prob- 
lem for such models. We note that in the present model the 
protoplanetary/protostellar companions, required as 'gates' 
would not be detectable spectroscopically in currently out- 
bursting systems, since the central star is swamped by disc 
emission during FU Ori outbursts. The detection of suit- 
able protoplanetary/protostellar companions (which induce 
radial velocity shifts of the central star of less than lkm/s) 
in T Tauri stars may however be within the reach of a new 
generation of high precision spectrometers. 

It is not clear at present what fraction of FU Orionis 
outbursts are of the rapid rise variety. Two out of the six 
confirmed FU Orionis systems are in this category, but this 
statistic needs to be interpreted with caution given the ob- 
servational bias towards identifying systems whose bright- 
ness has changed dramatically (Bell etal. 1995). Although 
this diversity of light curve behaviour may alternatively re- 
sult from intrinsic differences in disc properties, it is clearly 
undesirable to have to postulate variations in such funda- 
mental properties as the viscosity parameter, for example, 
from object to object. 
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In the present picture, as in the scenario proposed by 
Bell & Lin (1994), the epoch of FU Orionis outbursts is ter- 
minated when the mass input to the outer disc falls below 
the critical value at which the thermal ionisation instability 
can be triggered somewhere in the disc. As the mass supply 
dwindles, ta also falls (equation 10) so that the satellite is 
pushed to progressively smaller radii: depending on its de- 
tailed equation of state and structure it may or may not be 
tidally disrupted before finally being accreted by the central 
object. (The issue of dynamical tidal instability however re- 
quires further exploration (analogous, for example to that of 
Hjellming & Webbink 1987 in the case of composite stars) 
for objects with the complex equations of state of protoplan- 
etary/protostellar companions). If, on the other hand, the 
disc surface density is abruptly reduced (due, for example, 
to the clearing action of stellar winds) then such a satellite 
may avoid being accreted, being instead stranded at its in- 
stantaneous location when the disc was dispersed. Detailed 
calculations of tidally coupled thermally unstable discs are 
required in order to ascertain whether satellites are likely 
to be stranded in this way, and thus whether such objects 
- similar to that recently discovered in 51 Pegasus (Mayor 
and Queloz 1995) - should be commonly found in Weak Line 
T Tauri stars and post- T Tauri stars. 
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